Mass dispersions from giant dipole resonances using 
the Balian-Veneroni variational approach 



J. M. A. Broomfield and P. D. Stevenson 

Dept. of Physics, University of Surrey, Guildford, Surrey, GU2 7XH, United Kingdom 

Abstract. The Balian-Veneroni variational approach has been implemented using 
a 3-dimensional time-dependent Hartree-Fock (TDHF) code with realistic Skyrme 
interactions and used to investigate the mass dispersions from giant dipole resonances 
in 32 S and 132 Sn decaying through particle emission. The fluctuations obtained are 
shown to be quantitatively larger than the standard TDHF results. 
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1. Introduction 

The time-dependent Hartree-Fock (TDHF) approach can be used to determine the 
expectation values of single-particle observables, such as fragment mass, in nuclear 
reactions and decays but is known to underestimate the fluctuations in these values 
[1]. This is due to the 1-body nature of TDHF and the fact that it neglects 2-body 
correlations. This problem has previously been studied by Balian and Veneroni [2, 3, 4], 
who derived a general variational theory for the determination of expectation values, 
correlations and fluctuations. They found that, given the state of a system described, 
at the time to, by the 1-body density matrix, p (to) (a Slater determinant satisfying 
p 2 = p), the fluctuation, AQ, in a 1-body observable, Q, at some later time t±, is given 
by 

(AQ BV f\ ti = lim ^Tr [p (t ) - a (t , e)} , (1) 
where a (t, e) is a 1-body density matrix related to p (t) through the boundary condition 
cr (t u e) = exp (ieQ) p (h) exp (-ieQ) , (2) 

and where the time evolution of p (t) and a (t, e) is given by the usual TDHF equation. 
This result is significantly different from the standard TDHF result 



(AQ TDHF )\ = {(Q 2 )-(Q) 



l2\ /^i\2 



fi 



= Tr[Qp{t 1 )Q(l-p{t 1 ))], (3) 

in that it depends on the initial time, to, with the final time, t\, entering only through the 
boundary condition (2). The other key feature of this result is that it contains, through 
(2), the observable Q such that this method is specifically tuned to the determination 
of the fluctuation of the observable of interest. 

A practical implementation of (1) requires that a Hartree-Fock calculation be 
performed to determine the initial state, p (to). The system is then excited by a suitable 
external excitation, and a TDHF calculation performed from to — * t\ to determine 
p(t\). This is used to obtain a(t\,e) using (2) and a second TDHF calculation is 
then performed with the TDHF code run backwards, t\ — > to, to obtain a(to,e). The 
transformation (2) and the second TDHF calculation must be repeated for a range of 
values of e to allow AQbv to be determined by extrapolation to e — > 0. 

The large number of computations required to evaluate (1) and the complexity of 
these calculations means that only a handful of calculations have been performed using 
this method and those calculations which have been performed have used simplified 
interactions and made use of symmetries (either spherical [5], or axial [6, 7]) to render 
the problems tractable. However, modern advances in computing power mean that 
this approach can now be implemented using fully 3-dimensional TDHF codes with full 
Skyrme interactions [8, 9, 10, 11]. 

We consider the mass dispersion in a bounded region of space around a giant 
dipole resonance (GDR) which decays through particle emission and calculate the mass 
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(number of nucleons) in the nucleus according to 

N(R C )= £ Uf\K{f)\ 2 0{R c -\f\), (4) 

m<e F J 

where R c is the cutoff radius used to define the bounded region of space. 

The nucleus was excited by multiplying the ground state wavefunctions from the 
HF calculation by a dipole boost given by 

B D (x, y, z) = exp {iFC (A x x + A y y + A z z)) (5) 

with 

C = */A 1 (6) 

V 47r i _|_ exp (^fx 2 + y 2 + z' 2 ) 

and where, for protons, F — 1/Z, and for neutrons, F — — 1/(A — Z), where A is the 
atomic mass number of the nucleus under investigation and Z is its charge. A x , A y and 
A z determine the strength of the boost applied to the nucleus. 

Written in terms of the single particle wavefunctions (1) becomes [6] 

(AAW)\ = A-hm^, (7) 



/(e)= £ /df|^ ro (t o ,f,e)|0 n (to,f)>| 



2 



m,n<ep 



The wavefunctions |0 n (£)) were obtained from the results of a static Hartree-Fock 
calculation, whilst the wavefunctions \ip m (t,f,e)) result from the backwards TDHF 
calculations and are related to the wavefunctions \<f) n (t,f)) through the boundary 
condition 

V> (ti, f, e) = exp (ie0 (R c - |f|)) (t u f) . (8) 



2. GDR in 32 S 

We consider first a GDR in 32 S calculated using the Skyrme interaction with the SLy6 
[12] parametrisation. All calculations were performed in a cubic model space of size 
32 x 32 x 32 fm discretised in steps of 1 fm. The initial HF calculation gave a 32 S ground 
state with a total binding energy of 260.36 MeV (compared with the experimental value 
of 271.78 MeV [13]) and a prolate deformation = 0.11). 

At the beginning of the dynamic calculation the ground state wavefunctions were 
boosted in accordance with (5) and with A x = A y = A z = 112.5 fm -1 . The simulation 
was allowed to run from an initial time to — fm/c to t% = 250 fm/c in steps of 0.2 
fm/c. The emitted nucleons were reflected back from the boundary of the box and 
would, were the simulation allowed to run long enough, re-enter the region occupied by 
the de-exciting nucleus causing unphysical interactions. An analysis of the density and 
of (N (R c )) as a function of time was used to verify that the number of nuclei in the 
nucleus had stabilised well in advance of the time t\ and that the radiated flux had not 
had enough time to be reflected back and to interact with the nucleus. 
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Figure 1. (a) The dipole moments (Q x , Qy and Q z ) plotted as a function of time 
for a GDR in 32 S. The difference between Q x and Q y and Q z is consistent with a 
calculation using a prolate deformed ground state where x is the long axis, (b) ANg V 
plotted as a function of e and extrapolated back to s = 0. The standard TDHF result 
(calculated at t\ and independent of e) is shown for reference. 



The dipole moments, Q x , Q y and Q z , were obtained as a function of time using [10] 

Q, = ( (xf ) _ {x n )) j (9) 

where i = 1,2,3 denotes x, y and z and (xf ) and (xf) are the expectation values for 
position calculated using the proton and neutron single particle states respectively. This 
is shown in figure 1(a). Due to the prolate deformation of the 32 S nucleus, the Q y and 
Q z values are identical and differ from the Q x values. The periodicity of Q x , Q y and Q z 
allow the excitation energies for the oscillations along each of the three primary axes 
to be estimated. In this instance we obtain, for Q x , a period of ~ 71 fm/c giving an 
excitation energy E x « 17.5 MeV and, for Q y and Q z , a period of ~ 68 fm/c giving an 
excitation energy E y « E z m 18.3 MeV. 

The final state gave (N) = 26.65 with AN^ DHF = 4.08 using R c = 8 fm which 
represents the emission of ~ 5 nucleons. R c was chosen so that the bounded region fully 
enclosed the nucleus but omitted, as much as possible, the extended (or dissipated) 
components of the wavefunctions. The transformation (8) was then applied and the 
TDHF code was run in reverse. This process was repeated for e values in the range 
0.05 < e < 0.95 in steps of 0.05. At the end of each time-reversed calculation the 
fluctuation, AiV^y (e), was estimated using (7). These values were plotted (see figure 
1(b)) and a straight line was fitted to the linear section of the graph and extended back 
to e = to obtain AN^ V = 5.92 which represents a 20% increase in AiV using the 
BV approach compared with the standard TDHF result. This graph is typical of those 
obtained using this approach and is linear for larger values of e increasing asymptotically 
as e — > due to the 1/e 2 term in (1). Often, as in this case, the curve decreases for 
intermediate values of e where the reduced value of e means that the transformation 
(8) only has a small effect making the numerator in (1) numerically approximately zero 
and dominant over the e 2 denominator. 
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Figure 2. (a) The dipole moments (Q x , Qy and Q z ) plotted as a function of time for 
a GDR in 132 Sn. Q x and Q y and Q z as a result of the ground state being spherical. 
The shoulder at w 40 fm/c is a consequence of the 8 fm cutoff radius, (b) ANg V 
plotted as a function of e and extrapolated back to e = 0. The standard TDHF result 
(calculated at t\ and independent of e) is shown for reference. 



3. GDR in 132 Sn 

These calculations have been repeated for the doubly magic nucleus 132 Sn. All the 
calculations were carried out using the same model space and interaction as the 32 S 
calculation. The HF calculation produced an undeformed ground state with a binding 
energy of 1099.71 MeV (compared with the accepted value of 1102.85 MeV [13]). The 
ground state single particle wavefunctions were boosted at the start of the TDHF 
calculation in accordance with (5) and with A x = A y = A z = 600 fm -1 and the 
calculation was run from to — fm/c to t\ — 250 fm/c as in the previous calculation. 
The dipole moments were plotted as a function of time and are shown in figure 2(a). 
The graph shows Q x , Q y and Q z to be identical as expected for a spherical nucleus 
and gives the periodicity of the dipole moments as ~ 88 fm/c which corresponds to a 
resonance energy of ~ 14.1 MeV. This is close to the experimentally measured value of 
16.1 (7) MeV [14]. 

The standard THDF calculation gave, at the time t±, (N) = 121.02 and AN^ DHF = 
8.46 representing the emission of 11 nucleons. A series of transformations and time- 
reversed TDHF calculations were carried out as previously. The resulting graph, and 
linear fit, are shown in figure 2(b) which gives ANg V = 13.30 and represents a 25% 
increase in AiV compared with the standard TDHF result. 

4. Conclusions 

The Balian-Veneroni approach has been implemented for the first time using a 3- 
dimensional TDHF code with the full Skyrme interaction. Calculations have been 
performed for GDRs in 32 S and 132 Sn and have demonstrated that the BV approach 
does produce quantitatively larger results for the fluctuations of 1-body operators. This 
approach is now being applied to heavy ion collisions. 
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